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Abstract. We study the entanglement properties of a higher-integer-spin Affleck-Kennedy- 
Lieb-Tasaki model with quantum group symmetry in the periodic boundary condition. We 
exactly calculate the finite size correction terms of the entanglement entropies from the 
double scaling limit. We also evaluate the geometric entanglement, which serves as another 
measure for entanglement. We find the geometric entanglement reaches its maximum at the 
isotropic point, and decreases with the increase of the anisotropy. This behavior is similar 
to that of the entanglement entropies. 
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1 Introduction 

Quantum entanglement is a fundamental feature in quantum mechanics, and is a primary re- 
source in quantum communication and quantum computation [6, 12, 23, 37]. Entanglement has 
become an important tool to characterize quantum many-body systems (see [2] for example for 
a review). In one dimensional spin systems, typical quantifications of quantum entanglement 
are the Renyi entropy Sr(L,£) and von Neumann entropy S v n(L,£) of a subsystem A with i 
sites and environment B with L — £ sites (see Fig. 1) 

= logTr(p(L,l))^ SvN (L,£) = Hm S R (L,£). 
I — a a-tl 

Here the reduced density matrix p{L, i) is obtained from the density matrix of a ground state \^S>) 
by tracing out all spin degrees of freedom in the environment B 

The entanglement spectrum, i.e. the set of the eigenvalues of the reduced density matrix, deter- 
mines the entanglement entropies. For one-dimensional gapless spin chains, the generic behavior 
of the entanglement entropies has been analyzed [7] by use of the conformal field theory. The 
entanglement entropies scale logarithmically with the size of the subsystem, the prefactor essen- 
tially given by the central charge of the corresponding conformal field theory. 

On the other hand, gapful chains have been analyzed by investigating particular models. One 
of the most important models is the Affleck-Kennedy-Lieb-Tasaki (AKLT) model [1] which was 
introduced to understand the massive behavior of integer spin chains [13, 14]. The entanglement 
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Figure 1. Schematic pictures of the entanglement entropies (left) and the geometric entanglement (right) 
for the g-AKLT model with S = 3. 

entropies of the isotropic AKLT models have been investigated by examing the exact valence- 
bond-solid (VBS) ground state [9, 16, 17, 18, 22, 28, 31, 40]. For gapped systems which have 
finite correlation lengths, the entanglement entropies saturate at certain values when the size 
of the subsystems exceed certain lengths. The saturated values of higher rank and higher spin 
AKLT models are larger than the spin-1 AKLT model. 

Recently Santos et al. found surprisingly simple and useful formula for calculating the reduced 
density matrix for matrix product ground states [32, 33]. They applied it to the AKLT model 
of spin-1 and general integer spin S with quantum group symmetry (g-AKLT model) [3, 5, 
10, 19, 20, 24, 35], and another massive Kliimper-Schadschneider-Zittartz model [21] to study 
anisotropic effect. 

In this article, we study the entanglement properties of the q-AKLT model, following the 
results of [32, 33] and giving remarks and additional results. The more precise definition of the 
(/-AKLT model on an L-site chain with the periodic boundary condition is as follows 

25 

n=Yl E Cj(k,k + i)(Trj) k , k+1 , (2) 

fcez L j=5+i 

where Cj(k, k+1) > 0, and {ivj)k,k+Xi which acts on the fc-th and (/c+l)-th sites, is the U q {su{2)) 
projection operator from V$ (g) V$ to Vj, where Vj is the (2j + l)-dimensional highest weight 
representation of the quantum group U q (su(2)) [8, 15]. The valence-bond-solid (VBS) ground 
state of this hamiltonian T~L has a matrix product form [3, 24], which generalize the isotropic 
higher-integer-spin [4, 11, 36] and spin-1 g-deformed AKLT models [5, 19, 35]. We check that the 
entanglement spectra for t = 1 calculated from the formula of the reduced density matrix [32, 33] 
reproduce the one point functions originally derived in [3]. We achieve the finite size corrections 
of the entanglement entropies from the double scaling limit, which requires the second order 
term of the perturbation of the entanglement spectrum. We exactly calculate the finite size 
correction term of the von Neumann entanglement entropy S v n(£). 

Besides the entanglement entropies which characterize the bipartite entanglement, we also 
study the geometric entanglement, which is another kind of measure for entanglement, see Fig. 1. 
The geometric entanglement has been proposed as a measure for multipartite entanglement. It 
has been used to study quantum phase transitions [25, 26, 27, 28, 29, 30, 34, 38, 39], and 
has been measured experimentally recently [41]. Systems near criticality exhibit logarithmic 
divergences as the entanglement entropies. On the other hand, only a few analytic results are 
known for gapped systems. The geometric entanglement defined below can be regarded as the 
actual distance between the ground state of the system and the nearest fully separable state in 
the Hilbert space. 

We divide the L-site chain into N parties (L = Ni). Consider a pure quantum state of N 
parties |*} G H = ®f =l H^, where is the space of the ith party. The entanglement can be 
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quantified by maximizing the fidelity |A| between the quantum state \^f) and all the possible 
separable and normalized states of N parties |$) = (gi^l^M), G H®, 

($1*) 

|A max | = max |A|, A = -7==- (3) 

The logarithm of |A max | is taken 
£?(*) = - Log |A D 



v niitx| ) 



such that its value becomes zero when is separable or positive otherwise. The geometric 
entanglement per block is defined as the above quantity per party 

m = _ lim 

well defined in the thermodynamic limit. We evaluate the geometric entanglement for the spin S 
(/-deformed VBS state \^). We obtain the expression of the geometric entanglement for I — > oo 
and its finite size corrections with help of numerical calculations. For the evaluation of the 
entanglement entropies and the geometric entanglement, the spectral structure of the transfer 
matrix of the (7- VBS state in the matrix product representation [3, 24] will be helpful. 

This article is organized as follows. In Section 2, we briefly review the matrix product repre- 
sentation [3, 24] of the VBS ground state of the (/-AKLT model, which helps us for evaluating the 
entanglement entropies and the geometric entanglement. In Section 3, the finite-size correction 
terms of the entanglement entropies from the double scaling limit are calculated by perturba- 
tive analysis. We emphasize that the double scaling limits of the entanglement entropies and 
the leading term of the finite-size correction of the entanglement spectrum have been originally 
obtained by Santos et al. [32]. But we make Section 3 partially overlap their results so that 
this article can be self-contained and easy to read. In Section 4, we investigate the geometric 
entanglement with help of numerical calculations. Section 5 is devoted to the summary of this 
article. 



2 qr-VBS state 



In this section, we briefly review the matrix product representation of the higher-integer-spin q- 
VBS ground state and the spectral structure of the transfer matrix of the (/-AKLT model [3, 24]. 
We use the following notations. For a real number c we define its q analogue as 

M f ~ (f 

c = rr- 

q-q 

We also define the (/-shifted factorial and the (/-shifted binomial for n G Z>o as 

[n]! 



n 



i=l 

1, n = 0, 



[n]\[n-k]V A: -°'---' n > 
0, otherwise. 



The (/-VBS state [3, 24], which is the exact ground state of the (/-AKLT model (2), is expressed 
in the following matrix product form 



I*) = Tr[0i*52*"-*5£-i*5£]> 
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where g k is an (S + 1) x (S + 1) vector-valued matrix acting on the k-th site whose element is 
given by 



5fc (a,&) = (-l) s -V a+ ^ 5)(5+1)/ ^ 

=: h ab \S;b - a) k , (0<a,b<S). 



[S -a + b]\[S + a-b}\ \S;b-a) 



y\Oi) X y^\P)y Z 



Y. 



for vector-valued matrices 



The symbol * denotes the product A * B = 
A = {\a) xy } xya ndB = {\(3) xy } xy . 

We define gl by replacing each ket vector in the matrix by its corresponding bra vector: 

gl(a,b) = h ab k (S;b-a\. 

Let us set an (S + l) 2 dimensional vector space as 

W= C\ab)), 

0<a,b<S 

where {\ab}} | a, b = 0, . . . , S} is an orthonormal basis. We define an (S + 1) 2 x (5 + 1) 2 matrix G 
acting on the space W as 

G = g ] ®g, {{ab\G\cd)) = g\a, c)g(b, d) = 5 a - Ctb - d h ac h bd , 

which plays the role of a transfer matrix. 

In [3], the spectral structure of the G matrix was clarified, i.e. the eigenvalues of G are given 

as 



A n = (-1)"([5]!) 2 



'25 + 1 
S — n 



n = 0,1, ...,S, 



with the degree of degeneracy 2n + 1, and thus the squared norm of the ground state is given as 
<*|*)= TrG L = ( 2n + 1 ) X n- (4) 

0<n<5 

The matrix G has the following block diagonal structure since {(ab\G\cd)) = for a — b ^ c — d: 

min(S,S-j) 

G= G ij \ G {j) eEndWj, Wj = C\i,i + j)). 

S<j<S i=max(0,-j) 

The size of each block is (S -\j\ + l)x(S- \j\ + 1). Each element of G® is 

((a, a + j\G (j) \c, c + j)) = q (*+<*3-S)(S+l) [ 5 _ a + c]! [5 + a _ c ]i 





"5" 




S 




~S~ 




S 




a 




a + j_ 




c 




c + j 



We construct intertwiners among the 2S + 1 blocks G^ (j = —S,...,S). This helps us to 
construct eigenvectors of each block from another block with a smaller size. 
Let us define a family of linear operators {Ij}i<\j\<s as 



Ij : Wj -> 11} 1 (j > 0), 
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{a,a + j - l\Ij\c,c + j) 



a [a + j] [S-g-j + 1] 

1 V bW-j + i] 



c = a, 



_ a _ j I [a] [S-a + 1] 



U) [s-j + i] 



c = a — 1, 

0, otherwise, 
Wj -> Wj+i (j < 0), {{a -j -I, a\Ij\c - j, c)) = {{a, a-j- l\I-j\c, c - j)). 



By direct calculation, one finds that the matrix Ij enjoys the intertwining relation IjG^ = 
GV-Vlj (1 < j < S), IjG® = GV +r >Ij {-S < j < -1). Each block has a simple 
(nondegenerated) spectrum 

SpecG w = {\e}\j\<£<s, 
and the corresponding eigenvectors are given by 



|A„», 



0<i<S-n 



{n+1)l [S-n]\[i + n]\[S-i]\ 
V [S]\[n]\[S-i-n)\[z\\ 



+ n)), n = j > 0, 



5]![n]![5-i-ra]![i]! 



0<i<5-n 

• • • ^«|A n )) n , 

/i-l^--2---^n|A„))_ n , 

The ^th-power of the G matrix is formally expanded as 

^ = (B E ,//\ U V\ . l A ™>)j j(( A n|• 



1 < j + 1 < n < 5, 
1< -j + 1 < n < S. 



3 Finite size correction of the entanglement entropies 

In this section, we examine the finite-size correction of the entanglement entropies by studying 
the reduced density matrix. Recently, the following simple formula for the reduced density 
matrix (1) was found [33] 



p(L,£) 



K(L - £)K(£) 



TV G L 

where the U K matrix" is defined as 

K(£) = MG e 
with a linear map M 



(6) 



M(\ab))((cd\) 



ac 



(7) 



The reduced density matrix (6) is an (S + l) 2 x (S + l) 2 matrix, from which the rank of the 
density matrix is equal to or smaller than (S + l) 2 . We study the reduced density matrix by 
combining (6) and the spectral structure of the transfer matrix G reviewed in the last section. 
Here we introduce some notations and make some general remarks. We define 



K, 



1 



-A*(|A„»;;«An|) 
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so that the K matrix and the reduced density matrix are written as 

K(£) = >iKn, 
0<n<S 

^ fan + lU* S X n>^nK n >K n . (8) 

t-' J< ' n 0<n<S 

0<n<S o<n'<s 

One observes that K n , K{£) and p(L,£) enjoy the same block diagonal structure as G: 

Kn= K%\ K(£) = K^(£), p(L,£)= p^(L,£), 

-s<j<s -s<j<s -s<j<s 

K® , {£) , (L,£) 6 End Wj , 

since {(ab\K n \cd)) = for a — b ^ c — d. Note that M (7) does not always map a matrix acting on 
a sector Wj to a matrix acting on the same sector. The spectrum of p(L, £) is, of course, given by 
the union of the spectra of p^(L,£)'s. Due to the symmetry {{ab\p(L,£)\cd)) = ((ba\p(L,£)\dc)), 
we have the degeneracy 

Spec p^(L,£) = Spec (L,£). 

3.1 Double scaling limit 

We first review the double scaling limit [32, 33] 

p = lim lim p(L,£), p® = lim lim p®(L,£). 

£—>oo L—>oo ' I— >oo L— >oo 

Noting the form (8) and |A n /Ao| < 1 (n = 1, . . . , S), we find the reduced density matrix becomes 
diagonal 

q 2(a+b~S) 

p = K K , ((ab\p\cd)) = 5 ac 5 bd (=: 5 ac 5 bd p ab ). (9) 

[S + 1\ 2 

The entanglement spectrum is, of course, given by the diagonal elements of p, i.e. {p a b|a, b = 
0, 1, . . . , S}. 1 We notice that the degree of the degeneracy of the eigenvalue rg+ip 1S S — \k\ + 1. 
For example, the spectrum for S = 2 is given as 

1 



Specp (2) : p 2 



[3] 2 ' 



q (i) 1 r 

Specp v > : pox = -oiwfo, P12 = T^fo, 
rPr [3] 2 

c (o) 1 1 "7 4 

Spec^J: Poo= (? 4^p, Pii = P22 = ^p, 

Q (-1) 1 9 2 

Specp^ ^ : pxo = ^p, P21 = pjp 

Specp(- 2 ): P20 = t4- ( 10 ) 



[3] 



One can calculate 



0<a<S 
0<i><S 



lr This notation is different from that in [321 
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Figure 2. The von Neumann entropy in the double scaling limit 5 v n (12)- 
Then we achieve the entanglement entropies in the double scaling limit [32, 33] 
Sr 



Log P = Log 



1 — a 



I- a b [a][S + l] a '' 



a i g —l a S+l a_ g -(S+l)\ 

S M = 2Logq[^^ I -(S+l) q s+1 _ q +2Log[5 + l], 



q-q- 



(12) 



see Fig. 2 for the von Neumann entropy in the double scaling limit. In particular, when q = 1, 
the spectrum is totally degenerated 



Pab 



(13) 



(S + l) 2 ' 

and the entropies become 

S R = S v N = 2Log(S+l), 

which agree with the case of the open boundary condition [16, 22, 40]. On the other hand, in 
the limit q — > 0, only one eigenvalue survives poo = 1, Pab = (a + 6 > 0), and the entropies 
become zero. 



3.2 Finite-size correction 

We examine the finite-size correction of the entanglement entropies. We first take the limit 
L — > oo 



p{£):= lim p(L,£) = K V K n K £ n , 

L— >oo z — ' 

0<n<S 

L— ¥oo 



p^(£):= lim pV){L,e)=K® £ i^M, 



(14) 
(15) 



0<n<S 

j^, and then consider the case £ = 1 and the behavior of the entropies for 



with K n = and tnen consider tne case £ = l and tne benavior ot tne entropies tor l — > oo. 
Fig. 3 provides plots of the spectrum Spec p{£) of the reduced density matrix (14), i.e. the union 
of the spectra Specp^)(^)'s of (15), and the von Neumann entropy for S = 2 with q = 4/5. 
For £ = 1 and L — > oo the reduced density matrix becomes 



p(l) = K (MG)/X . 
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9 10 11 12 




3 4 5 6 7 



Figure 3. The entanglement spectrum of p{€) (left) and the von Neumann entropy (right) for S = 2 
with q — 4/5. The lines are drawn for I £ K with replacement n l n — » \n n \ and — > Hn^n^ 1 in (14) 
or (15). In the left figure, the dashed, dotted and solid lines correspond to Spec p^ ±2 \£) , Specp^ 1 ^^) 
and Spec p^(i), respectively. 



The eigenvalues of p(l) become zero except 25 + 1 ones, which is pointed out for S = 1 and 2 
in [32, 33] 



$a+b,k 
0, 



[S+\k\]\[S 



[25 + 1]! 



S-\k\ 

E' 

c=0 



(S+2)(2c+\k\-S) 



S) 


S 




S 




c 




c+\k\ 



a x b = 0, 



otherwise. 



(16) 



For example, the spectrum {p a b(l)}ab of p(l) for 5 = 2 is given as 

1 

M' 

i + <? 8 



Specp (2) (l) : 


P02(l) = 


Specp (1) (l) : 


Poi(l) = 


Specp(°>(l) : 


Poo(l) = 


Specp (_1) (l) : 


Pio(l) = 


Specp^ 2 )(l) : 


P20(l) = 



q 2 (l + q*)[b] 



6 „10 i „12 



1 - q z + 2<f - <? iU + 
9 4 (l + 9 4 )[5] 
l + o 8 



g 2 (l + 9 4 )[5]' 
l 

[5]" 



Pia(l) = 0, 

Pll(l)=0, P22(l)=0, 
P2l(l) =0, 



Let us show (16). We consider the submatrix p( k '(l), k > 0. The case for < is similar. 
By direct calculation, we find the matrix elements of (S — k + 1) x (5 — k + 1) submatrix p^ (1) 
are given by 



«C C + fc|p« d)|c + j, C + j + *» = (_l)i,2c +fe -5 +(2c+ i +fe -S )( 5 + l ) t^lJ^- 







5 




S 




S 


c 




c + j_ 




c + k 




c + j + 



(17) 



The rank of /^^(l) is 1, since the element of p^ k \l) (17) has a form yl c x B c+ j. Thus, only one 
eigenvalue of p^ k \l) is nonzero, which is given by 

S-k 



Tr (fc) p( fe )(l) = ^ ((c, c + fc|p«(l)|c, c + fc» 



c=0 
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Figure 4. The von Neumann entropy of the one-site subsystem 5 v n(1)- 



[S + k]\[S-k]l 
[25 + 1]! 



S-k 
c=0 



(S+2)(2c+k-S) 



S) 


~S~ 




5 




c 




c + k 



(18) 



from the fact that the other eigenvalues are all 0. The expression (18) is actually identical to 
the one point functions derived in [3]. In particular, when q = 1, the non-zero eigenvalues are 



degenerated as p a &(l) 



i 



25+1 



(a x b = 0), and we have 5 R (1) = 5 vN (l) = Log(25 + 1) [32]. One 



observes the monotonicity of the von Neumann entropy 5 v n(1) while < q < 1, see Fig. 4. 

We turn to the behavior of entropies for I — > oo. Noting again the form (8) and \k u \ < 1 
(n = 1, . . . , 5), we find 



p(e) = P + 4K Ki + o(4), 



oo. 



We denote the eigenvalue of p{£) by p a b(£) corresponding to p a b (9) when £ — > oo. Since the 
density matrix p in the double scaling limit is a diagonal matrix, it is not difficult to perform 
perturbative calculation. Noting \n\\ 2 > \k<2\ > \k^\ > • • • , we find 



Pab{t) = Pab + r a bK{ + t ab nf + o(«f ) , £ } OO, 

{{cd\ K Q K X I aft)) «o6| K Q Ki \ cc 



(19) 



r a fe = ((ab\K Ki\ab)), 



<<ab 



£ 

(c,d)#(a,6) 
0<c, d<S 



Inserting (5) into r a b and t a {, defined above, we have 
qd+b-S g a + fe - 35 [3](i _ q 2a 

r a b 



2a+2 _|_ q2S+2\M 



-,2b 



? 2fe+2 + q 2S+2^ 



tab 



[5 + 1] 

/ q 2(a+b-S) [3][2] ' 

1 [5 + 2][5 + l][5] 



(^-l)2 [5 ][5 + l][5 + 2] 



(20) 



q~ 4 [5 - a + 1] [a] [5 - 6 + 1] [6] - [5 - a] [a + 1] [S - b] [b + 1] 



(21) 



The first-order term (20) has been originally obtained in [32] (see equation (59) of [32] by 
changing the indices p = 5/2 — a, v = 5/2 — 6 and redefining q —> q 1 ^ 2 ), where the characteristic 
length is given by £ = Log([s+2]/[s]) ' ^ e a ^ so ca l ciua t e d the second-order term (21) which is 
needed for seeing the finite-size correction of the von Neumann entropy. 

For example, the spectrum {p a b(ty}ab — > oo) for 5 = 2 (which is shifted from (10) as (19)) 
is given as 



Spec^ 2 ^) : 



1 

w 



[2] 
[3] [4] 



"l! 
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g2[ 3 ]2 g 2 [3][4] i (l-g4)[ 4 p 

^ _ (l-g 2 )[2] f g 2 [2] 2 M 

[3]2 [3] [4] 1 (W)[ 4 p i' 

SpecpW(^) : + + g 2 (1 _ 9 4 )[4] 2 K i £ ' 

J_ (l-g 2 ) 2 [2] < _ [2] n g 4 g 2 [2] , g 6 [2] 2 « 

[3] 2+ 9 2 [3][4] Kl [4] Kl ' [3] 2 + [3][4] Kl (1 - g^]*" 1 ' 

bPeCP {l >- g2 [3] 2+ g 2 [3 ] [4 ] K l+ (l_ g 4 )[4] 2 K l » 

g 2 (l-</ 2 )[2] , g 2 [2] 2 



/' - - > ,-21 



Specp ( - 2) (£) 



[3]2 [3] [4] 1 (l-g*)[4f» 
1 [2] .j 



M 2 [3] [4] 11 

where we omit the symbol +o(k 2 ^). 

The Renyi entropy is expressed by p a b, r a b and i a fc up to the order of k 2 ^ as 

SR(^) = -^-LogTr(p(i)) B = -^-Log £ (p a ^)) f 
1 — a 1 — a _ z — L 

1 



0<a<S 
0<b<S 



Lo g {p(i + ^ + rKf)+ ( K f)} 



1 — a 

1 - a' 1 l-a\" 2 



5 R + :r~~~ /Cj + — - — (t — + o(«*) , . - x . i 22 ) 



where 



0<a<S 0<a<S \ 

0<6<S 0<6<S 

with P defined by (11). 

By tedious but straightforward calculation, one finds 

a[3] / [g + 2][a5]-[ a (S + 2)][5] \ 2 
[5] [5 + 2] V (9- T 1 )^ + !)][« + !] J ' 



T = a 



mm 



(q- q -i)2[S][S + 2] 



[2(q - 1)] /([(S + 2)(q + 1)][5] - [5 + 2][5(a + 1)]) 



[2] V [a + 2}[a + l][a(S + 1)] 



a-1 / [20S + 1)] 2 _ 2[2(5 + l)][a][(a + l)(g + l)] [2][a][(a + 2)(S + 1)P 2 
2 V[2][5 + l] 2 [S+l][a(S + l)][a + l] [(5 + l)a][a + 2] 

Then we find 

S vN (/) = S vN + - + , ^ °°, (23) 

where the coefficient of k\ vanishes. Since the leading order term is k 2 ^, the characteristic length 
is 2£. We find the coefficient of depends on the anisotropy parameter q but is independent 
of the spin value S. 
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As discussed in [32], the perturbation fails for the isotropic case due to the degeneracy (13), 
but the entanglement spectrum can be written by linear combinations of K n 's and has the same 
spectral structure for the transfer matrix G. For example, for S = 2, we have 



Spe C/ 9 (2) 0?) : 
Spec/9 (1) (^) : 
Specp (0) (^) : 
Specp< -1 >(i) : 



— + — 
2 6 



— + — 
2 6 



— + — 
2 6 



— + — 
2 6 



— + — 
2 6 



K 



3 + 2 



1 

3 + 2 

1 k\ 

- + — 
3 2 



6 

•jm*, PC Ms-fa. 



6 

5/ej 
~6~ 



1 A 



» + Ki H 

3 V 3 1 3 



where no higher order term is needed. In [32] the finite-size corrections of the entanglement 
entropies for q = 1 were calculated as 



S R (£) = 2 Log(5 + 1) - -anf + o{nf) , 
S yN (£) = 2 Log(S + 1) - + o( K f ) , 



oo, 



oo, 



which agree with the limits q — > 1 of (22) and (23). 



4 Geometric entanglement 

In this section, we evaluate the geometric entanglement, which is another kind of measure of 
entanglement. We divide the chain into A?" parties (L = N£), and each of the N parties to be 
contiguous blocks of i spins S. When N is large enough, the following expression for the fidelity 
|A max | (3) has been shown for PT-symmetric matrix product ground states \^f) in [25, 26, 29] 

\r1\ 2N 

i A -i 2 =i^Mw (24) 

where |d| 2 is the quantity 

\d\ 2 = max (Aux|G^| Aux), |Aux)= ^ x a x* b \ab)). (25) 

{Xi}: £ |X,P = 1 0<a<S 

i=0 



0<b<S 



Performing the maximization (25), one obtains the fidelity |A max | which finally leads to the 
analytic expression for the geometric entanglement £(£) 

= _ Log |A max | 2 ^ 

For convenience we set x% = rie^~^ 6i (r, > 0, flj 6 R), and write x* = r*e^~^ di if the setting 
{xi = x'} maximizes (Aux|G^| Aux). 
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0.0 0.2 0.4 0.6 0.8 1.0 0.0 0.2 0.4 0.6 0.8 1.0 



Figure 5. The geometric entanglements for £ = 1 (left) and 2 (right). In each graph, the line corresponds 
to S = 1 (26), and the markers x (5 = 2), □ (S = 3) and + (S — 4) are plotted based on numerical 
calculations. 



4.1 Spin-1 

We calculate the geometric entanglement for 5 = 1. By direct calculation, we have 

( AuK|rf| Aux) = + W< + (f + 2(-l/ - WW 

[3]" , rH-l)^ I- ,4 



+ 1 l2T + ^2^ ,|a;i 
Inserting r 2 = 1 — Tq, we get 

( Aux|^| Aux) = ^L{W ~ (-DO (d - 1)4 + 2 q rl) + |» + ^ 
where #j's do not appear. Thus we find 

0<g<l: rS = l, \d\ 2 = + ^f^, 



1 + g 2 1 + g 2 



2 _ 3* , 



g = l: |d| 2 = ^- + 



2 2 



*>1= ^ = 0, |dP = i^ + L^, (26) 

where \d\ 2 = ( Aux |G £ | Aux) is independent of {xi} at the isotropic point q = 1 [25], and the 
choice of r* changes discontinuously at this point. (We will see that this kind of "degeneracy" 
occurs for the higher spin case.) Inserting these forms and (^1^) = [3]^ + 3(— 1) L into (24), we 
finally achieve the geometric entanglement 

< q < 1 : £(£) = Log (l + q 2 ) - Log (l + g 2 (-[3])- £ ), 

q = l: £^) = Log2-Log(l + (-[3])^), 

g > 1 : £(£) = Log (1 + q 2 ) - Log (q 2 + (-[3])^), 

which generalizes [25]. The entanglement entropy takes its maximum at the isotropic point, 
decreases with the decrease of the anisotropy parameter q and finally becomes £{£) = at 
q = 0, see Fig. 5. This behavior of the geometric entanglement is similar to the entanglement 
entropies. In the limit I — > oo, we have 

£ = Log(l + q 2 ). 
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4.2 Spin-2 

Let us consider first the isotropic case, where we have 
(Aux|G £ | Aux) 

- 1 \t + 1 \i.L 1 \t (1 ~ (rg - r 2 ) 2 ) 2 - 4r r 2 r 2 (cos(fl - 26 1 + fl 2 ) + 1) , , ^ 
" 3 + 2 1 + 6 2 2 ^ " A ^ 

(^i ~ 2r r 2 ) 2 - 4r r 2 r 2 (cos(fl - 2^ + fl 2 ) - 1) e . 

with A = 40, Ai = -20 and A 2 = 4. When t is odd (resp. even), X\ < A 2 (resp. X[ > A|). 
Using the first (resp. second) form, we find 

£odd: r«-r« = 0, cos(^ - 20\ + 9' 2 ) = -1, \d\ 2 = h l + h e 2 , 



even: rf - 2r V 2 " = 0, co S (0S - 2B\ + 0' 2 ) = 1, \d\ 2 = \x £ + \x[ + \x l 2 . 

6 2b 



In the anisotropic case, thanks to the form 

+ 

we have 



A 2 2 

(Aux \G'\ Aux) = ^ ^ ri 4 r2 (A^ - X e 2 ) cos{d - 2d x + 2 ) + (independent of {0*}), 



£odd: cos(0*-2^ + ^) = -l, 
£even: cos(6>S - 26>J + 6>*) = 1, 

which is the same as for the isotropic case. We use help of numerical calculations (see Fig. 6 for 
1 = 1 and 2), which indicates that 

t odd : r* = 5 i0 for q < 3 q*; r* > 0, r* = 0, r* > for q* < q < 1, 

£ even : r* = 5{q always maximizes (Aux \ G \ Aux). 

One observes that the geometric entanglement with I odd is not completely monotonic while 
< q < 1, see Fig. 5. The set {r^r*} is obtained by 



'0 
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in the case where £ is odd and q* < q < 1 



. -(1 ~ g 2 )(l + Q 4 )K + 2[3]Aj - (1 + 3g 2 + + 
'" ' W (l-g 2 ) 2 (l + g 4 )A^ + 4 (7 2[3]A £ 1 - (1 + 2^ + 6^ + 2^ + 58)^- 



The transition point q* is obtained by solving (27) = 1, which approaches 1 as I — > 00. The 
set {r*} for > 1 is obtained by replacing ro ri and (/ — > 1/g. Under the assumption r* = <5jo, 
we have 



12 



£(£) = - Log Jit = - Log(Aux Aux)| ri=(Ji( 



2 , 4N t /, , g 2 (l + g a + ? 4 M , « 



= Log(l + ^ + <f) - Log ^1 + » - - - - + 1 + q]J . 
for < (7 < 1 and sufficiently large I. 
4.3 Spin-3 

For the isotropic case, we have 

( Aux|^| Aux) = ^ + + \M + 1a| - y(A? - Al), 



.4 



(rir 2 - 3r r 3 ) 2 + l{r\ - \F?>r\r?) + 2(r 2 - v / 3r r 2 ) 2 

- 2rir2{3ror3(cos(a + /?) — 1) + 2v / 3?"ori(cosa — l) + 2v / 3r 2 r 3 (cos/3-l)} 



where a = 9 - 26 l + # 2 and p = e l - 26 2 + 9 3 . When £ is even, X{ > A|. Thus ( Aux |6?'| Aux) 
is maximized by 

rf - \fh\r\ = rf - V3r^ = 0, cos a* = cos /3* = 1, 

When I is odd, the candidates of {0{} that maximize ( Aux|G^| Aux) for given {rj} are 

3r 2 r§ — 4r 2 r| — 3rgrf 3rpr 2 — 4r 2 r| — 3r 2 r§ 

COS OL — ^= - , cos p = ^= - , 

4v3ror^r 2 AySrir^r?, 



cos (a + l3) = - ' 1 g _ 1,2 7 ' 2 ' 3 , (28 



3 r 2 r 2 _ 4 r 2 r 2 _j_ 3 r 2 r 2 

6r rir|r 3 

or cosa = — 1, cos/3 = — 1, cos(a + /3) = l, (29) 

or cosa = l, cos/3 = — 1, cos(a + /?) = — 1, (30) 

or cosa = — 1, cos/3 = 1, cos(a + /?) = — 1. (31) 

Inserting (28), we get 



9 '« 2 1 n 2 1 2 ^ 



^ = 4 -( 3r + 2r I +r 2 9 



and thus we find 



2 



3r « 2 + 2rf + rf = *-, |d| 2 = + -/ 2 + (32) 
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0.0 0.2 



0.6 0.8 



Figure 7. Numerical calculations of the sets {r*} which maximize |<i| 2 for S — 3 and 
(right). The markers Q, A, x and □ correspond to r', r', r* and r*, respectively. 



1 (left) and 2 



We end up achieving the same value \d\ 2 (32) for (29)-(31). For example, inserting (29), we get 

2 



A = \~ ( 3r o + 2r l + r l ~ I ) - (vWi - n r 2 



+ V3r 2 r 3 ) . 



The maximization for the anisotropic case with i odd is more complicated than S = 2, see 
the numerical result in Fig. 7. We expect that 

i odd : there exist q°, q" , q"' and q"" such that 

r* = 5 i0 , < q < q', 

r* = r* = 0, r' > 0, r* > 0, q* < q < q", 
r* > 0, i = 0,1,2,3, q"<q<q 9 ", 
r\ = 8 a , q"' < q < q"", 

[r' =r' 2 =0, r\ > 0, > 0, 9 "" < 9 < 1, 
£ even : r* = <5^o always maximizes (Aux \ G e \ Aux). 

We also expect that these transition points q*,...,q**** approach 1 as £ — > oo. Under the 
assumption r' = 5io, we have 

\d\ 2 



£{£) = -Log\^ = - Log(Aux|G £ | Aux)| ri=5i0 



= Lo M -'[4]-Lo g (l + ?aM + -^ 
for < q < 1 and sufficiently large t. 



+ 



14 £ 
9 «3 



g 2 + 9 4 (l-g 2 + g 4 )[5] 



4.4 General case 

We consider the maximization of (Aux|G^| Aux) for general S. As we observed in the previous 
subsections, we expect that, for given q < 1, 

i odd: there exists I* such that the set {r* = 5io} maximizes (Aux \G l \ Aux) when £ > £*, 
I even: the set {r* = 5jo} always maximizes (Aux|G^| Aux). 
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£ 1.0 




Figure 8. The geometric entanglement in the limit I — >• oo. 



Since the term 



E <f r a 
a=0 



dominates in (Aux \ G | Aux) for I — > oo, we have r* — > Sio, 



o«Ao|Ao» 

which supports the above assumption. Inserting {r^ = 5^}, we have 



\d\ 2 = (Aux|G £ | Aux) 



n=Sio 



E 

k=0 



o((Afc|Afc)) ' 



and find 



£(£) 



Log 



' s 



^ o((Afe|Afe)) 



\k=0 



Log(q S [S+l}) 



[5 + 2] 



Ki + o Ik 



i 



oo. 



(33) 



Here we used the norm (4) of the (/-deformed VBS state |\&) and the norm of the eigenvectors 
of the transfer matrix [3]. In the limit i — > oo, we have the geometric entanglement £ = 
Log(g' s [5 + 1]), which takes the maximum Log(S' + 1) at q = 1 and approaches as q — > 0, see 
Fig. 8. The monotonic behavior while < q < 1 is similar to the entanglement entropies. 

The isotropic point is a special case where the choice r{ = 5oi or rj = 5si does not always 
maximize \d\ 2 for S > 2 even if I is large, as we saw for S = 2 and 5 = 3. Thus the asymptotic 
form (33) is no longer valid at the isotropic point. 



5 Summary and discussion 

In this article, we studied some entanglement properties of the higher spin g-AKLT model with 
the periodic boundary condition from the matrix product representation of the g-VBS ground 
state. We exactly calculated the finite-size correction terms of the entanglement entropies 
by the perturbative calculation for the spectrum of the reduced density matrix. We found 
that the first-order correction term of the Renyi entropy vanishes by taking the limit a — )■ 1. 
This requires the second-order perturbation of the entanglement spectrum for calculation of the 
finite-size correction of the von Neumann entropy. It would be interesting to extend the study 
of entanglement properties to various generalizations, the entanglement entropies with multiple 
blocks (see [31] for the isotropic spin-1 case), for example. We also investigated the geometric 
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entanglement. The geometric entanglement in the limit I — > oo decreases with the decrease of 
the anisotropy parameter q while < q < 1. This property is the same as the entanglement 
entropies. Under an assumption which is based on numerical results, we calculated the finite-size 
correction of the geometric entanglement. 
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